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ABSTRACT 

We study Schwinger pair creation of charged particles due to the inhomoge- 
neous electric field created by the thin electron layer at the surface of quark stars 
(the electrosphere). As suggested earlier, due to the low photon emissivity of 
the quark-gluon plasma and of the electrosphere, electron-positron pair emission 
could be the main observational signature of quark stars. To obtain the electron- 
positron pair creation rate we use the tunnelling approach. Explicit expressions 
for the fermion creation rate per unit time per unit volume are derived, which 
generalize the classical Schwinger result. The finite size effects in pair production, 
due to the presence of a boundary (the surface of the quark star), are also consid- 
ered in the framework of a simple approach. It is shown that the boundary effects 
induce large quantitative and qualitative deviations of the particle production 
rate from what one deduces with the Schwinger formula and its generalization 
for the electric field of the electrosphere. The electron-positron pair emissivity 
and flux of the electrosphere of quark stars due to pair creation is considered, 
and the magnitude of the boundary effects for this parameters is estimated. Due 
to the inhomogeneity of the electric field distribution in the electrosphere and of 
the presence of the boundary effects, at high temperatures (T > T cr « 0.1 MeV) 
we find a lower electron-positron flux as previously estimated. The numerical 
value of the critical temperature T cr depends on the surface potential of the star. 
We briefly consider the effect of the magnetic field on the pair creation process 
and show that the magnetic field can enhance drastically the pair creation rate. 

Subject headings: quark stars: electrosphere: electron-positron pair production 
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1. 



Introduction 



The Schwinger mechanism of pair production (Schwinger 1951), first proposed to study 
the production of electron-positron pairs in a strong and uniform electric field, has been 
applied to many problems in contemporary physics. Strong electromagnetic fields lead to 
two physically important phenomena: pair production and vacuum polarization. A strong 
electric field makes the quantum electrodynamic vacuum (QED) unstable, and consequently 
it decays by emitting a significant number of boson or fermion pairs (Greiner et al. 1985). 

For a spin 1/2 particle Schwinger 's predicted production rate per unit time and volume 
w is given by (Schwinger 1951; Soffel et al. 1982) 



where m e and e are the electron mass and charge, respectively, h is the transverse momentum 
and E is the (constant) electric field. In the present paper we use units so that H — c — 
k B = 1. In these units, e is equal to a l/2 and 1 MeV= 5.064 x 10~ 3 fm _1 . 

In most of the physical applications the proper time method introduced by Schwinger 
(1951) has been used to calculate the pair production rate. The real part of the effective 
action leads to vacuum polarization and the imaginary part to pair production. Though that 
method is conceptually well defined and technically rigorous, it is generally difficult to apply 
it to concrete physical problems such as inhomogeneous electromagnetic fields. Schwinger's 
result was generalized to electric fields E 3 = E 3 (x±), which depend upon either light cone 
coordinates x± = x 3 ± xq, but not upon both, in Tomaras et al. (2000). The form of the 
result is exactly the same as in the original formula given by Eq. (1). The case of electric 
fields depending on both x + and x_, E 3 = E 3 (x+,x-) was considered in Avan at al. (2003). 

An alternative approach to pair creation was initiated by Casher et al. (1979, 1980), 
who re-derived Schwinger's pair production rate by semi-classical tunnelling calculations. 
The boson and fermion pair production rate by strong static uniform or inhomogeneous 
electric fields was derived, in terms of instanton tunnelling through potential barriers in the 
space-dependent gauge, in Kim and Page (2002). 

The Schwinger mechanism for particle production in a strong and uniform electric field 
for an infinite system was generalized to the case were the strong field is confined between 
two plates separated by a finite distance by Wang and Wong (1988). The production rates, 
obtained by solving the Klein-Gordon and Dirac equations in a linear vector potential, can 
be expressed in an exact analytical form. The numerical evaluations of the production rates 
have shown large deviations from the Schwinger formula, thus indicating a large finite size 
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effect in particle production. The explicit finite size corrections to the Schwinger formula 
for the rate of pair production for uniform electric fields confined to a bounded region were 
calculated, by using the Balian-Bloch expansion of the Green functions, by Martin and 
Vautherin (1988, 1989). 

One important astrophysical situation in which electron-positron pair creation could 
play an extremely important role is the case of the electrosphere of the quark stars (Usov 
1998a,b, 2001; Page and Usov 2002). Quark stars could be formed as a result of a hadron- 
quark phase transition at high densities and/or temperatures (Itoh 1970; Bodmer 1971; 
Witten 1984). If the hypothesis of the quark matter is true, then some neutron stars could 
actually be strange stars, built entirely of strange matter (Alcock et al. 1986; Haensel et al. 
1986). For a review of strange star properties, see Cheng et al. (1998). 

Several mechanisms have been proposed for the formation of quark stars. Quark stars 
are expected to form during the collapse of the core of a massive star after the supernova 
explosion, as a result of a first or second order phase transition, resulting in deconfined 
quark matter (Dai et al. 1995). The proto-neutron star core, or the neutron star core, 
is a favorable environment for the conversion of ordinary matter into strange quark matter 
(Cheng et al. 1998). Another possibility is that some neutron stars in low-mass X-ray binaries 
can accrete sufficient mass to undergo a phase transition to become strange stars (Cheng 
and Dai 1996). This mechanism has also been proposed as a source of radiation emission 
for cosmological 7-ray bursts (Cheng and Dai 1996). Some basic properties of strange stars 
like mass, radius, collapse and nucleation of quark matter in neutron stars cores have been 
also studied (Glendenning et al. 1995, 1997; Cheng and Harko 2000; Harko and Cheng 2000, 
2002; Harko et al. 2004). 

The structure of a realistic strange star is very complicated, but its basic properties can 
be described as follows (Alcock et al. 1986). Beta-equilibrated strange quark - star matter 
consists of an approximately equal mixture of up u, down d and strange s quarks, with a 
slight deficit of the latter. The Fermi gas of 3A quarks constitutes a single color-singlet 
baryon with baryon number A. This structure of the quarks leads to a net positive charge 
inside the star. Since stars in their lowest energy state are supposed to be charge neutral, 
electrons must balance the net positive quark charge in strange matter stars. The electrons, 
being bounded to the quark matter by the electromagnetic interaction and not by the strong 
force, are able to move freely across the quark surface, but clearly cannot move to infinity 
because of the electrostatic attraction of quarks. For hot stars the electron distribution 
could extend up to ~ 10 3 fm above the quark surface (Kettner et al. 1995; Cheng and Harko 
2003; Usov et al. 2005). The electron distribution at the surface of the quark star is called 
electrosphere. The effect of the color-flavor locked phase of strange matter on the electric 
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field in the electrosphere was discussed by Usov (2004). 

Photon emissivity is the basic parameter for determining macroscopic properties of 
stellar type objects. Alcock et al. (1986) have shown that, because of very high plasma 
frequency near the strange matter edge, photon emissivity of strange matter is very low. For 
temperatures T « E p , where E p m 23 MeV is the characteristic transverse plasmon cutoff 
energy, the equilibrium photon emissivity of strange matter is negligibly small, as compared 
to the blackbody one. The spectrum of equilibrium photons is very hard, with uj > 20 MeV 
(Chmaj et al. 1991). 

The bremsstrahlung emissivity of the quark matter and of the electrosphere have been 
estimated recently in Cheng and Harko (2003); Jaikumar et al. (2004); Harko and Cheng 
(2005). By taking into account the effect of the interference of amplitudes of nearby in- 
teractions in a dense media (the Landau-Pomeranchuk-Migdal effect) and the absorption 
of the radiation in the external electron layer, the emissivity of the quark matter could be 
six orders of magnitude lower than the equilibrium blackbody radiation (Cheng and Harko 
2003). The presence of the electric field strongly influences the radiation spectrum emitted 
by the electrosphere and strongly modifies the radiation suppression pattern (Harko and 
Cheng 2005). All the radiation properties of the electrons in the electrosphere essentially 
depend on the value of the electric potential at the quark star surface. 

The Coulomb barrier at the quark surface of a hot strange star may also be a powerful 
source of e + e~ pairs, which are created in the extremely strong electric field of the barrier. 
In the case of the energy loss via the production of e~e + pairs, the energy flux is given by 
F± = e±n, with e± = m e + T and 



and Ate is the thickness of the region with the strong electric field (Usov 1998a, 2001). 

At surface temperatures of around 10 11 K, the luminosity of the outflowing plasma may 
be of the order ~ 10 51 ergs -1 (Usov 1998a,b). Moreover, as shown by Page and Usov (2002), 
for about one day for normal quark matter and for up to a hundred years for superconducting 
quark matter, the thermal luminosity from the star surface, due to both photon emission and 
e + e~ pair production may be orders of magnitude higher than the Eddington limit. Hence, 
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electron-positron pair creation due to the electric field of the electrosphere could be one of 
the main observational signatures of quark stars. 

It is the purpose of the present paper to consider the Schwinger process of pair creation in 
the electrosphere of the quark stars, by systematically taking into account the main physical 
characteristics of the environment in which pair production takes place. The electric field in 
the electrosphere is not a constant, as assumed in the Schwinger model, but it is a rapidly 
decreasing function of the distance z from the quark star surface. Therefore, in order to 
realistically describe the pair production process one must consider electron-positron creation 
in an inhomogeneous electric field. To study the pair production process we adopt the 
tunnelling approach and the fermion production rate in the electric field of the electrosphere 
is derived. The production rate is a local quantity, it depends on the distance to the quark 
star's surface and it is a rapidly decreasing function of z. The emissivity and energy flux due 
to pair creation at the quark star surface is also considered. An important factor which can 
reduce significantly the pair production rate is the presence of a boundary (the quark star 
surface). For electric fields perpendicular to the boundary there is a significant reduction 
in the magnitude of the pair production rate, which is also associated with an important 
qualitative change in the process, which becomes a periodic function of the distance to the 
quark star surface. 

The present paper is organized as follows. In Section 2 we review the basic properties of 
the electrosphere of quark stars. The electron-positron rate production in the electric field 
of the electrosphere is obtained in Section 3. The boundary effects are considered in Section 
4. In Section 5 we calculate the electron-positron pair flux of the electrosphere and compare 
our results with those obtained by Usov (1998a,b, 2001). A brief summary of our results is 
given in Section 6. 

2. Structure of the electrosphere of strange stars 

In the electrosphere, electrons are held to the strange quark matter (SQM) surface by 
an extremely strong electric field. The thickness of the electrosphere is much smaller than 
the stellar radius R ~ 10 6 cm, and a plane-parallel approximation may be used to study its 
structure (Usov et al. 2005). In this approximation all values depend only on the coordinate 
z, where the axis z is perpendicular to the SQM surface (z = 0) and directed outward. To 
find the distributions of electrons and electric fields in the vicinity of the SQM surface, we 
use a simple Thomas-Fermi model considered by Alcock et al. (1986) and take into account 
the finite temperature effects as discussed by Kettner et al. (1995). 
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This model can be solved exactly, and all physical quantities of interest (chemical po- 
tential, electric field etc.) can be expressed in an exact analytical form (Cheng and Harko 
2003; Usov et al. 2005). The chemical equilibrium in the electrosphere implies that the 
electron chemical potential /i e satisfies the condition —V (oo) + /i e (oo) = —V + /j, e , where 
V is the electrostatic potential per unit charge and V (oo) and /j, e (oo) are the values of the 
electrostatic potential and of the electron's chemical potential at infinity, respectively. Since 
far outside the star both V (oo) and /j e (oo) tend to zero, it follows that fj, e — V (Alcock et 
al. 1986). 

The structure of the electrosphere and the corresponding radiation processes essentially 
depends on the value of V q , the electric charge density inside the quark star. When the 
temperature of the quark star core drops below 10 9 K, the strange matter becomes superfluid. 
At this temperature quarks can form colored Cooper pairs near the Fermi surface and become 
superconducting. From the BCS theory it follows that the critical temperature T c at which 
the transition to the superconducting state takes place is T c = A/1.76, where A is the pairing 
gap energy (Blaschke et al. 2000). An early estimation of A gave A ~ 0.1 — 1 MeV (Bailin 
and Love 1984), but some recent studies considering instanton-induced interactions between 
quarks estimated A ~ 100 MeV (Alford et al. 1998). 

Strange quark matter in the color-flavor locked (CFL) phase of QCD, which occurs for 
A ~ 100 MeV, could be rigorously electrically neutral, despite the unequal quark masses, 
even in the presence of the electron chemical potential (Alford et al. 1998). Hence, for the 
CFL state of quark matter V q = and no electrons are present inside or outside the quark 
star (Lugones and Horvath 2003). 

However, Page and Usov (2002) pointed out that for sufficiently large m s the low density 
regime is rather expected to be in the " 2-color-flavor Superconductor" phase in which only 
u and d quarks of two color are paired in single condensate while the ones of the third color, 
and s quarks of all three colors, are unpaired. In this phase, some electrons are still present. 
In other words, electrons may be absent in the core of strange stars but present, at least, 
near the surface where the density is lowest. Nevertheless, the presence of the CFL effect 
can reduce the electron density at the surface and hence it can also significantly reduce 
the intensity of the electric field and the electromagnetic emissivity of the electrons in the 
electrosphere. Therefore, in order to describe the radiation properties of the electrosphere 
we assume that V q ^ 0. 

The Poisson equation for the electrostatic potential V {z, T) generated by the finite 
temperature electron distribution reads (Alcock et al. 1986; Kettner et al. 1995) 

' 1 Sr = ^{(V*-V°) + n\V-V,)T*], Z <0, (4) 
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d 2 V 4a , , n N 

— = - { V» + WV),z>0, (5) 

where T is the temperature of the electron layer, which can be taken as a constant, since 
we assume the electrons are in thermodynamic equilibrium with the constant temperature 
quark matter. In Eqs. (4)- (5), z is the space coordinate measuring height above the quark 
surface, a is the fine structure constant and V q /3Tc 2 is the quark charge density inside the 
quark matter. The boundary conditions for Eqs. (4)- (5) are V — > V q as z — > — oo and V — > 
for z — > oo. In the case of the zero temperature electron distribution at the boundary z = 
we have the condition V(0) = (3/4)V q (Alcock et al. 1986). 

The general solution of Eq. (5) is given by (Cheng and Harko 2003) 

V{*,T) = ^ (6) 

where zo is a constant of integration. Its value can be obtained from the condition of the 
continuity of the potential across the star's surface, requiring V q (0,T) = V (0,T), where 
V q (z,T) is the value of the electrostatic potential in the region z < 0, described by Eq. (4). 
Therefore 



1/31 
Zo = —\ — — arcsmh 
2 V anT 



V2ttT 
K(0,T) 
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The number density distribution n e of the electrons at the quark star surface can be 
obtained from n e (z,T) = V 3 /3n 2 + VT 2 /3 (Kettner et al. 1995; Cheng and Harko 2003) 
and is given by 

= ^l + a ^l2^ + ^)] 
3 sinh 3 [2^fT(z + z )] 

In the limit of zero temperature, T — > we obtain V(z) = a /(z + b), where a = 
a/37t/2q; and b is an integration constant, b can be determined from the boundary condition 
V(0) = (3/4) V^, which gives b = (4a /3V^). Therefore, in this case we find for the electron 
particle number distribution the expression: n e (z) = (l/37T 2 )ag/ (z + b) 3 . 

In the absence of a crust of the quark star, the electron layer can extend to several 
thousands fermis outside the star's surface. 

The strength of the electric field E outside the quark star surface is given by 
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and can be expressed as 



1 3 sinh 2 [2yffT(z + z )] V ; 



The ratio of the electric field of the electrosphere of quark stars and of the critical electric 
field E cr = m 2 e /e is represented, as a function of the distance from the quark star surface 
and for different values of the temperature, in Fig. 1. 



3. Electron-positron pair production rates in the electrosphere 

Let us consider a quantized Dirac field coupled to a classical external electromagnetic 
field described by the potential A^. The probability P to remain in the ground state, i.e., the 
probability of emitting no pairs, is given by P = |(0 0)| 2 , where S is the S matrix defined 
as S = T exp \i J Lid A x\ , T is the time ordering operator and Lj is the interaction Lagrangian 
(Soffel et al. 1982). The probability P can also be written as P = exp [— J w(x)d 4 x] , where 
w (x) = 2lmL e ff(x), where L e ff(x) is the one-loop effective Lagrangian density, which 
includes all orders in the external field, but neglects self-interactions of the matter fields. 
The quantity w(x) can be interpreted as the pair production rate per unit time and unit 
volume at the space-time point x = (x , Xi,x 2 , x 3 ) (Greiner et al. 1985; Martin and Vautherin 
1989). 

An alternative description of the pair creation process can be obtained by assuming 
that the vacuum decays if there exists an ingoing antiparticle mode which is at the same 
time an outgoing particle mode. An in-vacuum is defined via the ingoing (anti)-particle 
basis, _a; |in, vac) = + a k |in, vac) = 0, whereas the out-vacuum is defined by the outgoing 
states ~ a\ |out,vac) = + a k |out,vac) = 0, where a are the particle creation and annihilation 
operators, respectively (Soffel et al. 1982). The number of the outgoing particles in the mode 
k created in the in-vacuum is N kk > = (vac, in| + a~l + a k i |vac, in) = \p kt \ 5 (k — k'), where 
the coefficients j3 k i are the elements of the single-particle S matrix. Due to the validity of 
the energy and momentum conservation laws, one can generally choose (3 k i = f3 k 8 k i and thus 
the number of created particles in the mode k is simply given by N kk , = \p k \ S (k — k'). The 
delta function can be re-expressed as a delta function over the frequencies of the associated 
modes, (3 k i = T k 5 (uji — uj k ). Then the rule {5 (uji — uj k )} 2 — > (l/2n) 5 (ui — uj k ) J dt may be 
used to calculate a continuous rate of creation of particle-antiparticle pairs (Soffel et al. 1982; 
Kim and Page 2002), 
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The physical meaning of the particle creation process is the following: a wave packet 
of negative frequencies incident on an electric field oriented along the z— direction will be 
partly reflected by the electric field and partly transmitted to z — > oo as a wave with 
positive frequencies. This process is nothing else than tunnelling (Soffel et al. 1982; Greiner 
et al. 1985; Wang and Wong 1988). However, in the following we restrict the tunnelling 
probability to the transmission probability through the potential barrier but exclude any 
non-zero transmission probability above the potential barrier. 

In the WKB approximation the transmission probability T can be approximately given 
as ITI 2 = exp (—2 f, \k z \ dz) = exp (—2a), where a — f, \k z \dz and k z is the longi- 

II *\ J barrier I / ^ \ n J barrier I ^ I *> o 

tudinal momentum of the electron (Soffel et al. 1982; Wang and Wong 1988; Kim and Page 
2002). 

For a fixed frequency iv, the momentum k z of an electron in the electrosphere of the 
quark stars is given by 

k z = ^ m l + kl-[uj-V(z)}\ (12) 

where kj_ = k\ + k 2 y (Soffel et al. 1982; Kim and Page 2002). The mean particle production 
number (JV) is given by summing over all modes (Greiner et al. 1985), 

IN) = ! -2a dtduJ dxdk x dydk y 

[ ' J 6 2n 2tt 2tt ' 1 ' 

In order to estimate the electron-positron pair production rate in the electrosphere, we 
have to find first the thickness of the classically forbidden zone, or, equivalently, the limits 
of integration z± in the transmission probability. They can be obtained as solutions of the 
equation k z (z) — and are given by 

z± = \ \ A^arcsinh [ ^7== t , 2 ] - *b- < 1 1 ! 

Therefore, by taking into account the form of the electrostatic potential in the electro- 
sphere we obtain for a the expression 



™ i+k ^-^-i^Sk^\ )iz - (15) 



With the help of the transformation r] — {l/ sj m\ + k\^j |c<j — \/2ixT / sinh 2y/an/3T (z + z ) 
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we obtain the following integral representation for a: 




3tt ra\ + k] 



a/1 — rfdrj 



(16) 



■ V -l 



2tt 2 T 2 _|_ ( y/ m l+ k j 



-1] — 1 



The variation of the transmission probability |T| 2 = exp {—2a) is represented, as a 
function of the parameter ( = \Jm 2 + fcj/o; and for different values of the temperature 



and particle energy ratio 2-k 2 T 2 /uo 2 in Fig. 2. For large values of the particle energy uo » 
\/m 2 + £r[, ( — > and the transmission probability is equal to 1, |T| 2 
probability also increases with the temperature of the electrosphere. 



1. The transmission 



In the process of particle production in the electrosphere of quark stars due to the 
tunnelling from the negative energy state, a positive energy particle (an electron) is created, 
leaving a hole in the negative energy continuum, which can be taken to be an anti-particle 
(a positron) moving in a direction opposite to the direction of the created particle. The 
created pair of particles is characterized by the energy (frequency) uo. For such a pair, 
one cannot, strictly speaking, specify a particular point as the location where the pair is 
produced (Wang and Wong 1988). One can only say that the pair of particles begins to 
emerge between the point z_ and z+. Nevertheless, one can associate the point z, which 
is the solution of the equation uo = V(z), as the location in the vicinity of which a pair 
of particles is produced. With this approximate association, the energy of the produced 
particle is then identified by an approximate location, and the energy interval also can be 
approximately related to a spatial interval at which the pair of particles is produced via the 
relation duo = (dV/dz) dz = eE (z, T) dz. 

Moreover, we shall introduce polar coordinates in the momentum space so that k x = 
k cos 6 and k y = ksin9. Therefore the electron-positron pair production rate ri±, giving the 
number of electron-positron pairs created per unit time and per unit volume by the electric 
field of the electrosphere at the surface of quark stars, (N) / At Ax Ay Az, can be written as 



ri± (z, T) 



An 2 



f 

Jo 



eE(z,T) x 



k exp 



r 37r" m\ + k 2 
'MV 2 (z,T) 



i a/W 



V(z,T) 



drj 



2tt 2 T 2 



rk 2 



V(z,T) 



where s describes the spin degrees of freedom of the produced particles (s 
and s = 2 for fermions). 



dk,(17) 



1 for bosons 



- 11 - 



Eq. (17) has a general formal structure which is very similar to that of Eq. (1), 
describing the electron-positron pair production in a constant electric field. In both cases 
the pair production rate is proportional to the intensity of the electric field, multiplied by an 
integral over the transverse momentum of the particles. However, since in the electrosphere 
E(z, T) — > for large z, the particle production rate also naturally tends to zero and outside 
the electrosphere or in the regions with low electric field the particle production ceases. The 
variation of the pair production rate in the electrosphere, as a function of the distance z is 
represented for different values of the temperature in Fig. 3. 



By introducing a new variable ( = \Jm1 + k 2 /V (z, T), the pair production rate can 
also be represented as 



h(z,T) = —eE(z,T)V 2 (z,T)x 



/oo 



Cexp 



9 j^_ 2 f 1 ^^(Crj-iy'drj 



2a I 27r 2 r 2 



d(. (18) 



In order to find an approximate representation of the pair production rate we power 
expand the integrand in Eq. (16). In the first order in 77 we obtain 

V 1 - 1 ! 2 _ w yAw 2 + k\ ^ 



u 2 



("27r 2 T 2 + J< 

Therefore in the first order of approximation a is given by 

^3~7r ml + k 2 



+ l^j r] + O (r] 2 ) . (19) 




(20) 



By fixing the electron energy so that uo = V(z, T) and taking into account the expression 
for the electric field in the electrosphere, given by Eq. (10), it follows that the electron- 
positron pair production rate can be written as 



2 „, mx r l «, /3tt m 2 + A; 2 



2 f 00 
n±(z,T) » — e£(z,T) jf exp 
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4tt 2 



ei? (z, T) I exp 



47T 3 



e 2 E 2 (z,T)exp 



— 71 



E r . 



m 2 e + k 2 
eE(z,T) 



kdk 



-7T 



E(z,T) 



(21) 



This is exactly the leading term in the Schwinger formula for the pair creation in a 
constant electric field E , but with the constant electric field substituted with the inhomoge- 
neous, z-dependent electric field in the electrosphere, E — > E(z,T). Therefore, in the first 
approximation, it is possible to study electron-positron pair production in arbitrary electric 
fields in the electrosphere by simply substituting the constant electric field in the Schwinger 
formula with the inhomogeneous electric field. 



In the second order of approximation we obtain 

I3tt5 m 2 + k 2 



a 6 



1 1 



f + 7T 2 T 2 



2^ 



1 3n (ml + k 2 ) 
+ 3V« 



X 



7T 2 T 2 



3 

+ t- 



4^ + tt 2 T 2 8^ + 7r 2 T2 y 



(22) 



giving 



h± (z,T) 



8tt 2 



e 2 E 2 (z, T) exp 



25tt 



48A (z,T) 



[3A(^,T)]- 1/2 [l-erf( X (^,T))], (23) 



where X(z,T) = eE[l/V 2 + {1/3tt)V 2 /E 2 + {2/3tt 2 )V 6 /E 4 }, X = v / X{z,T)/3irE cr /E + 
5v^F/4v^A (2, T) and erf (x) = (2/^/tt) Jf exp (-m 2 ) d«. 



4. Boundary effects in pair production in the electrosphere 

The electron-positron pair creation in the electrosphere takes place very close to the 
surface of the quark star, which represents the boundary of the system. Therefore electron- 
positron pair creation is localized to the half-space z > 0, and the pair creation takes place in 
electric field localized to a bounded region in the space. Boundary effects in electron-positron 
pair creation by electric fields confined in a finite region of space have been previously 
investigated, with the general result that finite-size effects induce large deviations of the 
production rate from what one deduces from the Schwinger formula (Wang and Wong 1988; 
Martin and Vautherin 1988, 1989). 

In analyzing the surface effects on pair production there are two very different situations, 
which lead to very different results. The surface effects essentially depend on the orientation 
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of the electric field with respect to the boundary (Martin and Vautherin 1988, 1989). In 
the case of an electric field parallel to the boundary, the pair production rate per unit 
time and unit volume in the region z > is w (z) = Y^=i w n [1 — ex P {—eEz 2 /nir)], where 
w n = (e 2 E 2 /4n 3 ) exp (— nnE cr /E) is the Schwinger production rate in constant field with 
infinite extension. As one can see for z = the boundary effects cancels the Schwinger volume 
term, and thus electron-positron pair creation cannot take place on a boundary parallel to 
the field (Martin and Vautherin 1988, 1989). Generally, there is a significant reduction of the 
Schwinger pair production rate near a surface which is parallel to the electric field. For an 
electric field of the order of E = 5 x E cr , the reduction in the production rate at a distance 
of z — 197.5 fm from the parallel boundary is, for the case of the dominant term with n — 1, 
[1 — exp (— eEz 2 /7i)] = 0.33, while for z = 987.362 fm the reduction rate is 0.999. For the 
case of an electric field E = 50 x E cr , we have a reduction rate of 0.98 at z = 197.5 fm and 
1 at z = 987.362 fm. 



However, the case of the electric field parallel to the boundary is not relevant for the 
quark star electrosphere, in which the electric field, which is oriented outward, is perpendicu- 
lar on the surface of the star (the boundary). In the case of the electric field perpendicular to 
the boundary one can find the correction terms by using a method based on approximating 
the Green functions in terms of classical paths. The first order correction to the effective 
Lagrangian for particle production is (Martin and Vautherin 1989) 



^(1) 

J eff 



w-s/ " 



eE coth (eEs) 



exp 



2 ■ E 2 

-im p s-\-ie — z coth 
2 



eEs 



(24) 



The first order correction W\(z) to the pair production rate is given by W\ (z) 
21mL,fff(z), and for m e m can be written in the form (Martin and Vautherin 1989) 



wi(z) = 



-^e 2 E 2 I $ (s) sin (^z 2 s ) ds, 



where 



(lnf±±)' 




(25) 



(26) 



The ratio of the total production rate w = wo + w\ and of the Schwinger production 
rate w = (e 2 E 2 /4ir 3 ) exp (—iiE cr /E) is represented, for different values of the electric field, 
in Fig. 4. 

For large z the function ImL^^(z) has the asymptotic behavior 



ImL^Jz) 



1 sin(fz 2 ) 
'Stt 26 hi(f^ 2 ) + 



1 



In (f z 2 ) 



(27) 
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while for small z it is given by 



lmL%(z) 



e 2 E 2 eE 2 
+ 



48tt 8tt 2 



eE 



c, 



(28) 



where c = 1/3 - [s$ (s) - 1/3] « -0.1 (Martin and Vautherin 1989). 



From these equations and from Fig. 4 it follows that when the boundary is perpendicular 
to the field the pair production rate exhibits oscillations as a function of the distance to the 
boundary. No such oscillations occur when the electric field is parallel to the boundary. For 
the discussion of the physical origin of this phenomenon and its possible implications see 
Martin and Vautherin (1989). 

As we have shown in the previous Section, in the first approximation we can substitute 
the inhomogeneous electric field in the results derived for the constant field case. We shall 
follow this approach in the study of the effect of the boundary (the quark star surface) on 
the pair production in the electrosphere. Hence we assume that the boundary effect in the 
pair production rate is given by 

1 f°° 
Wl ( z ,T)^-—e 2 E 2 (z,T) J $(s); 



sm 



eE(z,T) 2 ■ 
z s 



ds, 



(29) 



which is a straightforward generalization of the constant field case. 

Therefore the total rate production n±) (z,T) = h± (z,T) + w\(z,T) of the electron- 
positron pairs in the electrosphere, with the boundary contribution included, is given by 



n 



4tt 2 



eE(z,T) V 2 (z,T) x 



/ 



a cxp 



V(z,T) 

1 



-2 




1 a/I — i] 2 (arj — 1) 1 dr\ 

+ H - 1) 2 



4tt 2 



e 2 E 2 (z,T) 



$ (s) sin 



2tt 2 T 2 
V(z,T) 2 



da — 



ds. 



(30) 



The ratio of the total pair production rate in the electrosphere, with the boundary effects 
included, and the Schwinger pair production rate in a constant electric field is represented, for 
different values of the temperature, in Fig. 5. As expected, the boundary effects significantly 
reduce the creation of electron-positron pairs. 



5. Electron-positron pair flux of the electrosphere 



For an electron in the electrosphere sitting near the Fermi surface the energy is given 
locally by e = fi e (z,T) — V(z,T). Since from the boundary conditions it follows that 
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(j, e (z,T) — > and V(z,T) — > for z — > oo, the energy of the electrons satisfy the con- 
dition e = 0. From this result it immediately follows that de/dz = 0. Electron and positron 
pairs are created with an energy e± = m e + T, which is greater than the energy of the 
electrons in the electrosphere, e± > e. 

We define the electron-positron pair emissivity Q± of the electrosphere as the energy 
created per unit time and per unit volume by the electric field, given by the product of the 
number of pairs multiplied by the energy of each pair. The mathematical definition of the 
electron-positron pair emissivity in vacuum is 

Q±(z, T) = £± Af ^ Az = e±n± T) = (m e + T) n ± {z, T) . (31) 

Q±(z, T) is general a local quantity, depending on the distance to the quark star surface 
and on the temperature in the electrosphere. The variation of the electron-positron emissivity 
is represented, for different values of the temperature and for a fixed quark star surface 
potential, in Fig. 6. In the case of a constant electric field the electron-positron emissivity is 
Q<£> (T) = e±w = (m e + T) (e 2 /4vr 3 ) E 2 exp (-nE cr /E). For an electric field E = 50 x E cr 
we obtain Q { ± ] (T) = 1.28 x (m e + T), which for T = 3.5 MeV gives (T) = 5.12 MeV 5 . 

The electron-positron pair flux of the electrosphere in vacuum is defined as 
F° (T) = - r Q±(z, T)dz = - r e±n ± (z, T) dz = ^ + T) f°° n± (z, T) dz. (32) 

K Jo TT J 7T Jo 

In the case of a constant electric field E = 50 x E cr , and assuming that the electroshere 
extends to up to d = 500 fm, we obtain (T) = (m e + T) (e 2 /4vr 4 ) E 2 exp (-nE cr /E) d w 
1 x (m e + T) , which for T = 3.5 MeV gives F± (T) w 4 MeV 4 . The variation of the electron- 
positron flux is represented, as a function of the temperature and for different values of the 
quark star surface potential, in Fig. 7. 

The variation of the electron-positron emissivity in the presence of a boundary perpen- 
dicular to the electric field is represented, as a function of the distance z to the star's surface, 
in Fig. 8. The boundary of the star strongly suppresses the pair production and induces 
an important qualitative change in the behavior of the emissivity Q. The electron-positron 
flux in the electrosphere, by taking into account the boundary effects, is represented in Fig. 
9. There is a significant effect in the electron-positron flux due to the suppression of pair 
creation by the quark star's surface. 

The possibility of the thermal contribution to the pair creation process was discussed 
by Gies (1999). In the first loop approximation there is no thermal contribution to the 
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imaginary part of the effective Lagrangian and therefore to the particle production rate. 
This means that the pair production rate in strong electric fields is basically independent of 
temperature, and the temperature dependence of this effect in the electrosphere of the quark 
stars is indirect and is due to the temperature dependence of the electron number density 
and of the corresponding electric field. 

Finally, we compare the generalized Schwinger pair production rates and the corre- 
sponding flux with the expressions proposed by Usov (1998a,b), and given by Eq. (2). We 
adopt a typical model for the electrosphere with an electron number density n e of the order 
n e = 2 x 10~ 5 fm~ 3 m 154 MeV 3 . The chemical potential of the electrons, /i e , which is also 

1 /3 

equal to the Fermi energy of the electrons 6f, is given by \i e = (37r 2 n e ) ~ 16.58 MeV. The 
number density of the pairs created by the electric field, which in this model is taken to be 
equal to the density of electronic empty states with energies below the pair creation thresh- 
old at thermodynamical equilibrium, is given by An± sov ^ = (3T//i e ) exp (— 2m e /T) n e 
27.8Texp (—2m e /T). By assuming that the electron spectrum is thermalized due to electron- 
electron collisions, the electron-positron pair flux from the electrosphere can be written as 
F 07sov) = 01673 x T s x ( me + T) exp (-2m e /T) x J(f), where f = 1.59/T. For T = 0.1 
MeV, F^ sov) « 7.4 x 1CT 9 MeV 4 , f| Usov) « 0.53 MeV 4 for T = 10 MeV and f| Usov) « 3.45 
MeV 4 for T = 40 MeV. As one can see from Figs. 7 and 9, there is a very large difference 
between the thermalized electron flux and the "pure" Schwinger flux due to pair creation by 
the electric field. 

However, this comparison is not correct in the concrete physical framework of the quark 
star surface, since in the Usov (1998a,b, 2001) model the electron thermalization time is used 
to calculate the flux, while the Schwinger flux is estimated in terms of the pair creation time 
in vacuum, which is of the order of At w 1/ (eE) 1 ^ 2 (Nikishov 1970). In order to compare 
our results with the previous estimations of the electron-positron pair rate at the surface 
of strange stars, we have to systematically take into account the number of available (free) 
electron states at the surface of the star and the thermalization effects in both models. 

To compare the two mechanisms we start with the discussion of the created particle 
number densities in the two models. In the Schwinger mechanism in vacuum, the particle 
number density is An± = h(z,T)At = h(z,T) / (eE) 1 ^ 2 , which, in the approximation of 
the constant electric field, can be written as 




(33) 



For a fixed value of the electric potential in the electrosphere V q , which for simplicity 
we take as the electric potential at the surface of the quark star, the electric field is given by 
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E = y/4/37rV q ^V£/2 + tt 2 T 2 , and the electron number density is n e = V; 3 /3tt 2 + V q T 2 /3. 
The temperature dependence of An± is only through the temperature dependence of the 
electric field E, which is significant only for values of the temperature so that T > V q . 

However, the number of free electron states at the quark star surface is given by 
An (Usov) = 3Texp(-2m e /T) (V 2 /3tt 2 + T 2 /3) . From a physical point of view, An { ± sov) 
represents the number of available quantum states for pair creation (Usov 1998a,b, 2001). 
Moreover, these states are not necessarily filled, because of the finite value of the electric 
field. Hence, the number of created electron-positron pairs is also finite, and this number can 
be less than the number of available quantum states An± sov \ Therefore, electron-positron 
pair creation is possible once the condition 



Ani Usov) > An° ± , 



(34) 



is fulfilled. Generally, these condition will be satisfied for temperatures so that T > T cr , 
where the critical temperature can be obtained from the condition 



3T cr exp 



2m P 



V 2 



T 2 

cr 

3tt 2 ~3~ 



471-3 



E{T cr ) 
E 

J — J CT 



3/2 



exp 



E r 



-7T 



E(T cr ) 



(35) 



The value of T cr depends on the electrostatic properties of the quark star surface. The 
critical temperature is represented, as a function of the electrostatic potential V q at the quark 
star surface in Fig. 10. The critical temperature is increasing with increasing V q . Hence, 
when T < T cr , all quantum states are filled, whereas for T > T cr the number of available 
quantum states may exceed the number of electron-positron pairs created by the electric 
field. 

Therefore, by taking into account the number of free electron states we define the 
electron-positron pair number density An± as 



An 



± 



Ani Usov) , 



T<T„ 



An° ± ,T>T a 



(36) 



The variation of the particle number densities in the two models is presented in Fig. 11. 
For the chosen values of V q the condition An± sov ^ > An^_ holds for temperatures T ^ 0.1 
MeV for V q = 5 MeV and T > 0.4 MeV for V q = 15 MeV, respectively. It should be noted 
that despite the number of available free electron states increases very rapidly with the 
temperature, due to the temperature independence of the Schwinger process, the number of 
the created electron-positron pairs is basically determined by a single parameter, the quark 
star surface potential. 
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In the model of Usov (1998a, 2001) the thermalized electron-positron pair flux is given 

by 

4 th) = A7rR 2 Ar E An^ sm) t; h \ (37) 

where R is the radius of the star, Ar E the thickness of the emitting region and is 
the characteristic time of thermalization of the electrons. For the expression t^ 1 ps 
(3/2tt) (a/0r) {T 2 /e F ) J(f), with £ = 2 \focpK (ep/T), was assumed (Usov 1998a, 2001). 
The function J(£) is defined by Eq. (3). 

In order to obtain a more realistic description of the thermalized electron-positron flux 
of the electrosphere which takes into account the inhomogeneities of the electron and electric 
field distributions, we assume that, due to its dependence on the Fermi energy Bp ~ // e = 
V(z, T), the thermalization time is also a function of the distance z to the quark star surface. 

Therefore we define the thermalized electron-positron flux from the quark star's surface 
in the generalized emission model as 



/oo q poo 

An ± (z, T) t£ (z, T) dz » ^R 2 ^T 2 J An ± (z, T) 



J 



2 /g V(z,T) 
V 7T T 



V(z,T) dZ ' 
(38) 



The variations of the thermalized electron-positron fluxes in the Usov model and in the 
generalized electron-positron emission mechanism are presented in Fig. 12. 

In the limit of low temperatures the condition T/V q — > holds with a very good approx- 
imation. By assuming that the number density of the e + — e~ pairs can be approximated 
by Eq. (33), in which all the parameters are estimated near the quark star surface z ~ 0, 
the electron-positron thermalized flux from the electrosphere of the quark stars can be rep- 
resented in an approximate form as 



F| th) ^ A-nR 2 Ar E x V q 



83 3 / 4 (2tt) 1/4 \V q 




il7im^c W {-^^) nr;(^) (-) ,...}. (39) 

For T = the electron-positron flux from the quark star surface is zero. Generally, F± 
is a function of the electrostatic potential at the quark star surface V q and of the temperature 
T. 



The presence of a surface magnetic field H can also enhance the pair production rate 
(Nikishov 1970). The magnetic field increases the pair production rate by a factor 8 H = 
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11H/E coth (11H/E). If H » E, there will be a significant increase in the pair production 
rate. The electric field of the electrosphere could be as high as E = 40E crit ^120 MeV 2 . 
On the other hand the estimated magnetic fields at the surface of the quark stars could be 
of the order of H « 10 15 - 10 17 G « 20 - 2000 MeV 2 (1G = 1.953 x 10~ 14 MeV 2 ). 

Magnetic fields with such high values may be present in very young quark stars. As- 
suming equipartition of energy, the energy of the differential rotation can be converted into 
magnetic energy, so that IQ 2 (AQ/Q) ~ (Att/3) R 3 (H 2 /87r), where / is the moment of in- 
ertia of the star, R its radius and Q and AQ are the angular velocity and the variation of 
the angular velocity, respectively. Therefore the magnetic field of a young quark star can be 
approximated as H « 10 4 (Att/tt) 1/2 MeV 2 . By assuming that AQ/Q m 0.03, we can obtain 
values of the magnetic field as high as H Pzt 2000 MeV 2 . Of course magnetic fields of such 
strength are not stable, because they will be pushed to and through the surface by buoyant 
forces and then reconnect (Kluzniak and Ruderman 1998). For a magnetic field of the order 
of H w 2000 MeV 2 we have S H « 53. 

Therefore strong magnetic fields can significantly increase the electron-positron pair 
production rate, and, consequently the luminosity of the electrosphere of quark stars. 

6. Conclusions 

In the present paper we have re-considered the electron-positron pair emission from the 
electrosphere of quark stars, as originally proposed by Usov (1998a,b), by pointing out the 
important role the boundary effects and the inhomogeneity in the distribution of the electric 
field may play in the pair creation process. At zero temperature, there are no available 
free energy states in the electron plasma at the strange star's surface. Therefore, at low 
temperatures T < T cr pa 0.1 MeV (corresponding to a quark star surface electric potential of 
V q = 5 MeV), the pair production mechanism by the strong electric field of the electrosphere 
is severely limited by the quantum effects and the exclusion principle specific to the Fermi- 
Dirac statistics. At high temperatures T > T cr pa 0.1 MeV, the pair creation rate is controlled 
by the electric field E and not by the temperature, because such a process is essentially a 
quantum process. Once the number of available electron states becomes higher than the 
Schwinger pair production rate, electron-positron pairs can be freely created by the electric 
field at the surface of strange stars. This happens at a critical temperature T cr , which strongly 
depends on the electrostatic properties of the quark star surface. The critical temperature 
increases with the increase of the electrostatic potential V q . At high enough temperatures, the 
pair creation process is almost independent of the temperature and is controlled exclusively 
by the electric field. On the other hand, the actual thermalized pair creation rate, which, 
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from an astrophysical and observational point of view is the most relevant quantity, depends 
on the temperature through the thermalized time scale. The number density of pairs can 
be accurately evaluated by using the Schwinger formalism and by taking into account the 
inhomogeneities in the electric field distribution. However, the boundary effects also induce 
large quantitative and qualitative deviations of the particle production rate from what one 
deduces with the Schwinger formula and its generalization for the inhomogeneous electric 
field of the electrosphere. 

Due to all these effects, we estimate that at high temperatures the energy flux due 
to e~ — e + pairs production could be lower than in the initial proposal of Usov (1998a,b). 
However, this flux could still be the main observational signature of a quark star. On the 
other hand, the presence of a strong magnetic field at the quark star surface may significantly 
enhance the electron-positron flux. 

The possible astrophysical and observational implications of the direct pair production 
effect will be considered in a future publication. 
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Fig. 1. — Ratio of the electric field of the electrosphere E and of the critical electric field 
E cr = rn^/e, as a function of the distance z (fm), for different values of the temperature: 
T = 0.01 MeV (solid curve), T = 5 MeV (dotted curve), T = 10 MeV (dashed curve) and 
T = 15 MeV (long dashed curve). In all cases V q (0,T) = 15 MeV. 
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Fig. 2. — Transmission probability |T| 2 = exp (—2a) for electron-positron pair creation in 
the electrosphere of quark stars as a function of the parameter ( = ^m 2 e + k\/oo for different 
values of the temperature energy ratio: 2it 2 T 2 /uj 2 = 0.1 (solid curve), 2ir 2 T 2 / lu 2 = 1 (dotted 
curve), 2it 2 T 2 /uj 2 = 3 (dashed curve) and 2ir 2 T 2 /u 2 = 5 (long dashed curve). 
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Fig. 3. — Ratio of the electron-positron pair production rate n± in the electrosphere and of 
the Schwinger rate w for a constant electric field E = 50 x E cr as a function of the distance 
z (fm) for different values of the temperature: T = 0.5 MeV (solid curve), T = 1.5 MeV 
(dotted curve), T = 2.5 MeV (short dashed curve) and T = 3.5 MeV (long dashed curve). 
In all cases the surface electrostatic potential of the quark star is V q = 15 MeV. 
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Fig. 4. — Ratio of the pair production rate w = Wq + W\ for a constant electric field, with 
the boundary effects included, and the Schwinger rate w , as a function of the distance z 
(fm), for different values of the electric field: E = 5 x E cr (solid curve), E = 25 x E cr (dotted 
curve) and E = 50 x E cr MeV (dashed curve). 
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Fig. 5. — Ratio of the pair production rate ri± = n± + w\ in the electrosphere of the 
quark stars, with the boundary effects included, and the Schwinger rate of pair creation in 
a constant electric field Wq, as a function of the distance z (fm), for different values of the 
temperature: T = 0.5 MeV (solid curve), T = 1.5 MeV (dotted curve), T = 2.5 MeV ( short 
dashed curve) and T = 3.5 MeV (long dashed curve). The value of the constant electric field 
used to calculate the Schwinger rate is E = 50 x E cr . In all cases the surface electrostatic 
potential of the quark star is V q = 15 MeV. 
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Fig. 6. — Electron-positron pair emissivity of the electrosphere of the quark stars, as a 
function of the distance z (fm), for different values of the temperature: T = 0.5 MeV (solid 
curve), T = 1.5 MeV (dotted curve), T = 2.5 MeV ( short dashed curve) and T = 3.5 
MeV (long dashed curve). In all cases the surface electrostatic potential of the quark star is 
V a = 15 MeV. 
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Fig. 7. — Vacuum electron-positron pair flux F± of the elect rosphere of the quark stars, 
as a function of the temperature T (MeV), for different values of the surface electrostatic 
potential V q : V q = 8 MeV (solid curve), V q = 10 MeV (dotted curve), V q = 12 MeV ( short 
dashed curve) and V q = 14 MeV (long dashed curve). 
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Fig. 8. — Vacuum electron-positron pair emissivity Q± of the electrosphere of the quark 
stars, with the boundary effects included, as a function of the distance z (fm), for different 
values of the temperature: T = 0.5 MeV (solid curve), T = 1.5 MeV (dotted curve), T = 2.5 
MeV ( short dashed curve) and T = 3.5 MeV (long dashed curve). In all cases the surface 
electrostatic potential of the quark star is V q — 15 MeV. 
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Fig. 9. — Vacuum electron-positron pair flux F±^ of the electrosphere of the quark stars in 
the presence of the boundary effects, as a function of the temperature T (MeV), for different 
values of the surface electrostatic potential V q : V q — 8 MeV (solid curve), V q = 10 MeV 
(dotted curve), V q = 12 MeV ( short dashed curve) and V q = 14 MeV (long dashed curve). 
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Fig. 10. — Critical temperature T cr as a function of the surface electrostatic potential of the 
quark star. 



-33- 




T(MeV) 

Fig. 11. — Temperature dependence of the electron-positron number density An± sm ^ in 
the Usov mechanism (solid curve) and the particle number density An± in the Schwinger 
mechanism for different values of the surface electrostatic potential: V q = 5 MeV (dotted 
curve), V q = 10 MeV (short dashed curve) and V q = 15 MeV (long dashed curve). For the 
Usov mechanism we have assumed that the Fermi energy of the electrons is Ef = 20 MeV 
and the thickness of the electrosphere is Ar E = 1000 fm. 



-34- 



> 

| 0.008 

CM 

0.006 
5- 0.004 
^ 0.002 







0.1 0.2 0.3 0.4 

T (MeV) 



0.5 0.6 



Fig. 12. — Temperature dependence of the thermalized electron-positron flux F± h ^ /4irR 2 in 
the Usov mechanism (solid curve) and in the generalized electron-positron emission mech- 
anism for different quark star surface electrostatic potentials: V q — 5 MeV (dotted curve), 
V q = 10 MeV (short dashed curve) and V q = 15 MeV (long dashed curve). For the Usov 
mechanism we have assumed that the Fermi energy of the electrons is Bp = 20 MeV and the 
thickness of the electrosphere is Ate = 1000 fm. 



